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ABSTRACT 


It Is shotm that the flux fron lo fluctuates id.th an aaplltude of 

~10 percent on tine scales of every 28 seconds, 40 alnutes, night-to-night, 
and perhaps year-to-year. Such behavior Is found to be the result of random 
fluctuations for vhlch the mean square fluctuation varies Inversely with 
frequency for constant bandwidth measurement. It Is suggested that the theory 
developed for thermionic emission from barium oxide cathodes In vacuum tubes 
Is applicable to this situation. If this Is so, the fluctuations In the flux 
from lo's volcanoes may be caused by diffusion of hot convective cells onto 
the surface of lo. 

It is speculated that long-term fluctuations may furnish a means by 
which the lo volcanlsm can shut down and conserve energy. Thus the 
discrepancy that may exist between measurements of the current heat flow from 
lo and calculations of tidal dissipation may be resolved. 

Tests for rapid flickering at 10 )im showed no fluctuations greater than 1 
percent. This agrees td.th the prediction of the author's "flow model" theory 
(Slnton, 1982) In which the 10 ^m volcanic thermal emission arises from 
cooling of "old" flows. 

Running Head: lo Flickering 

Editorial Correspondence: Louise H. Good, Editor 

Institute for Astronomy 
2680 Woodlawn Drive 
Honolulu, HI 96822 
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I. INTRODUCTION 

The variability of Che 5 compoaenC on lo'a thennally emitted radiation 

has been demonateted In a number of papers (Wlttebome et al., 1979* 1982; 
Slnton, 1980; Morrison and Telesco, 1980; Slnton^el. » 1983). In the last 
named paper. It was demonstrated that at least some of the variability 
occurred in events or outbursts which have a duration of at least several 
hours. These occur somewhat at random, but seem to occur more frequently on 
the trailing hemisphere of lo and when lo la In the active sector of the 
Jovian magnetosphere. The cause of the variability Is attributed to 
variability of thermal emission from volcanic eruptions. 

In February 1981 while observing lo, 1 became aware that, according to 
the CRT display of the data acquisition system, 1 was spending fully as much 
time to achieve a specified S/N ratio on lo as I was on Ganymede and Calllsto, 

even though lo Is several times brighter at 4.8 pm than the other two 

satellites. It became evident that lo was “flickering" or changing brightness 
In a time comparable to the 28 seconds required for an Individual pair of 
integrations. Subsequently, data were taken In such a manner as to enhance 
the detection of rapid variations. After a determination of the spectrum of 
the rapid fluctuations, which showed that the spectrum was quite broad. It was 
realized that the previously reported nlght-to-nlght variability of lo and the 
rapid fluctuation might be part of the same phenomenon. The thesis of this 

paper is that the power spectrum of the fluctuations Is closely fit by an 

inverse proportionality to frequency. 

We discuss the nature of the rapid fluctuations In Section 11. After an 
Initial discussion of how the lo fluctuation can be separated from system 
noise, e statistical test of the reality of lo fluctuation la made, and this 
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Is followed by e determination of the power spectrum of the rapidly sampled 
data. A cursory reader can skip this section entirely and go to Section III, 
but Section II lays the foundation for the statistical analysis of the rapid 
fluctuation. 

In Section III, the fluctuation measurements at the various time scales 
are combined. It Is found that an assumed reciprocal frequency dependence of 
the fluctuation power Is In agreement with the observations. A physical otodel 
that predicts such a broad spectrum Is suggested In Section IV. In Section V 
it Is speculated that fluctuations of very long period may act to conserve the 
heat on lo and that the magnitude of the present volcanlsm Is atypical. A 
negative test for rapid flickering at 10 (im Is reported In Section VI. 


II. STUDIES OF THE RAPID 10 FLUCTUATION 


a. Determination of the Total Fluctuation 

The measurement of the amount of fluctuation of a variable source In the 
pres'ince of a nearly comparable measuring system noise Is not ai easy task. 
Only the total noise can be detected. To obtain an accurate value of the 
total noise a relatively large number of observations have to be made. The 
standard error estiniate of a signal measurement Is customarily determined 
from 


^ S 2 

1 h 


(n-l) 


1/2 


( 1 ) 


where Sj^ Is the 1th measurement of the signal and n Is the total number of 
data points. Normally one is primarily Interested In the signal measurement 
and is only secondarily Interested In the value of a to set bounds on the 
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value of S. If, however, one wants to compare the value of a for a 
fluctuating source with that for a noovarlabla source, then an accurate value 
of a Is required. The estimate of the standard error of a Is (Cameron, 1960) 


s s 2 s V 

1 ^1 n U ®1^ 


1/2 


( 2 ) 


2 (n-1) 

and the ratio of the standard error estimate to Its uncertainty Is 


2. - / 2 (n-1) 

<J 

a 


(3) 


Thus, In order to know a to 10 times Its uncertainty, 51 measurements have to 
be made. Normally, at the 2.2-m telescope an Integration requires 10^ and a 
pair of Integrations (the telescope is beam-switched by the data system 
between each) Is required for a signal determination. Settling time 
requirements lead to a total Intexrval of 28 seconds. Thus more than a 1000^ 
are required for a measurement of a good to 10 percent. 

The problem Is exacerbated In the presence of noise In the measuring 
system. Because the two sources of fluctuation are uncorrelated, the observed 
fluctuation Is given by 

- o. 2 + jy 2 (4) 

b 8 

and 

o ■ / (j2 - o. 2 (5) 

s b 


where Og Is a measure of the fluctuation of the source and ob 1 b the 
additional fluctuation of the system. The system noise can be determined by 
observing a star which Is assumed not to fluctuate. If Ob - then by the 
usual method of combining errors, the ratio of Og to Its uncertainty Is 
[ 2(n-l)/5]^/2lf n measurements are made of both the fluctuating source and the 
system noise. For a 10 percent uncertainty, n must be 251 for e^ch. 
Fortur.ately, In many cases the fluctuation of lo was substantially larger than 
the observed fluctuation of stars. 
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b« Determination of the Sytem Noise 

In the previous analysis we assumed only one source of noise In the 
telescope-detector system. Actually there are several, but for our purpose 
they may be treated as only twot one Is Independent of the signal level and 
one proportional to the signal level. The first may comprise detector noise 
and fluctuation of the telescope background and sky emission. Note that the 
number of photons from the sky and telescope are many orders of magnitude 
greater than those from the object. Thus the term that Is dependent on the 
square root of the signal Is Inslgnf leant. The signal-dependent noise sources 
are composed of guiding errors, seeing excursions, and scintillation caused by 
the atmosphere. At the signal level of lo, both the signal-dependent and 
signal-independent components are Important. To measure these and to 
determine the total system noise, 0 |), at lo's signal level, a number of stars 
spanning a wide range of brightness have been measured. We did not have 
available a single star having exactly the brightness of lo. For each 
observed star we have determined the nolse-to-slgnal ratio, N/S, for a single 
pair of 10® Integrations. This was generally determined from 16 pairs of 
integrations. The error of this quantity is, by (3), 1//30 of Its magnitude. 
For each observation, the total observed noise from the two components is 
given by the sum of their squares. This may be expressed as 

n 2 - Nb^ + C s2 (6) 

whert C^/2 is the proportionality constant for the signal (S) dependent term. 
Dividing by s2 gives us 

N 2 

n 2 ”b + C (7) 


It is convenient to express the signal In magnitudes on the right side of (7) 
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by putting S - 10'*^*'* M where M le the magnitude of the atar. Extracting the 
square root gives us the most useful form 



( 8 ) 


The constant D Is then We found C and D for each night by applying 

linear regression to the data In the form of (7) and where the Independent 
variable is 1/S^. However, the uncertainty of each N^/S^ varies as the 
brightness of the star and so a regression analysis that was weighted 
acrording to the S/N was used. Since similar results were obtained for the 
several nights of each observing run, all of the data In each run were 
combined. Figure 1 shows the data In the form of (8) with the least-squares 

I 

fitted curve. Error bars are not shown In these plots. With the exception of 
the discovery run, 22-25 February 1981, n was always 16 or an integer multiple 
of 16. Each data point plotted In panels b, c, and d of Fig. 1 represent 16 
measurements and so by Eq. (2) the error bars extend (1/30)^/^ of the value of 

the point. The curves are all excellent fits to the data over a range as much 

as eight magnitudes. There are, however, differences from one observing run 
to another. Between February and March 1981 the background noise Increased 
substantially, while better seeing apparently reduced the fluctuations of 
bright stars. A different detector system was Introduced in 1982. Because of 
its better quantum efficiency and the use of 6-arcsecond Instead of a 
9-arcsecond aperture, the noise at the faint star end in the 1982 data is 
substantially reduced. The noise at the bright star end is the same or 
reduced despite the smaller aperture. The reason for this Is apparently that 
the new system has a much more uniform response over the entrance aperture. 

The agreement at the faint end for the three runs In 1982 (see also Fig. 2) is 

most remarkable. Differences between runs are only a few percent. 
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c. Statistical Test of lo Fluctuation 

The measuremenCt of lo were made In exactly the same way as those of the 
stars. It Is clear from the location of the lo points on the plots of Fig. 1 
that a considerable number lie substantially above the curves from the star 
measurements. To see whether lo Is really flickering we have applied the 
F-test for the significance of the ratio of the larger variance (lo) to the 
smaller (stars). We used the critical F value for 15 degrees of freedom (n-1) 
for lo and ~400 degrees for the stars since typically the star fluctuation was 
determined from a combination of many observations. Using the double-tailed F 
critical value of 3.4 for 95 percent confidence, we have determined new curves 
in Figs. 1 and 2. If lo's N/S ratio lies above the appropriate curve, we have 
95 percent confidence that lo was flickering during that measurement. We can 
also view the data in Figs. 1 and 2 in a different light. There are a total 
of 8 lo values above the 95 percent level. By random chance (since there are 
91 lo determinations) we should expect less than 5. The stars appear to have 
approximately these numbers of high values. 

It is clear that many lo points lie on the least-square curve or within a 
standard deviation. Ue conclude that while lo sometimes flickers, it does not 
always do so, as evidenced by the plot in Fig. 2. 

Table I gives the part of the observed noise due to lo's fluctuation 
derived from (5). From the data in the table, we have formed the rms value. 
The rms is preferred over the straight average, not only because it is a more 
accurate appraisal of the fluctuations, but because theory, developed in a 
later section, requires it. In deriving the rms value, we have assumed that 
when the lo fluctuation is less or equal to that determined for stars, the 
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■quare root of the negative ouaber In (5) la laagloary and Its aquare reducea 
the Bum. The final reault of the lo raa flickering amplitude from the 91 aeta 
of measurementa of Table I la 0*060 ± 0*005 relative to Ita total algnal* 

d. Power Spectrum of the Rapid Fluctuation 

We subaequently aought to oeaaure the power apectrum of the flickering 
during one lun that waa available on the IRTF* Pig* 3 ahowa aeveral aelected 
time series of data In which we believe that lo waa flickering. Alao ahown In 
this figure are algnal meaaurementa of two atara whoae brightneaa la very 
similar to that of lo* Most of the IRTF measurementa were made at an 
integration of two seconds lo each of the two positions of the telescope 
(l*e>, five times as fast aa the 2*2-v data) since, at the time, we believed 
it was important to press toward the highest possible time resolution* As It 
turned out, most of the data do not have a fluctuation that la greater than 
the star data. Power spectra formed from the IRTF data of lo did not show any 
repeatable peaks In the spectrum. They may, however, show a trend toward an 
increase at lower frequencies, but It la not at all clear that this trend is 
not also present In the stellar data. 

We were more successful in determining the power spectrum from 2.2-m data 
we already had in hand and from some we subsequently obtained, albeit at 
longer sampling times. We selected 21 of the 91 data sets that have the 
highest F values so that we have a high degree of confidence that lo was 
flickering while these data were taken. The mean power spectrum of these data 
is shown as the solid dots lo Fig. 4. In computing this spectrum, we combined 
some contiguous sets of 16 or more data pairs. The computed power spectra 
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were "blDned" into 10 frequenclas of width 1.7 nHs bwglnalng at 1.7 and 
extending to 17 aHx. The naan apectrua whowa a narked decay toward higher 
frequenciea, and tha higheat flicker power ia found at tha loweat frequency. 

Becauae we aelectad data aeta in which lo had greater than avarage 
flickering, the power apectrun derived by the above proceaa ia expected to 
have a greater than average flicker power at each frequancy. To detemine the 
average apectrum, the valuea in the varioua frequency bine have been acaled by 
the ratio of the nean variance of the 21 aaoplea to the mean variance of the 
91 aamplea of Table I. It ia theao acalod valuea that are ahown in Pig. 4. 

In performing any Fourier apectrun analyaie one nuat alwaya aaaure 
oneaelf that the apectrun ia not aliaaed. The data aanplea that go into the 
spectrum of Fig. 4 are all iron differencing two 10-aecond integrationa apaced 
two seconds apart. Thus higher frequencies that would cauae aliasing are 
removed by the integration proceaa. Thu effective filtering corresponds to a 
sine function attenuation of offending frequencies. 

Since the peak of the flickaring appaara to be at atlll lower 
frequencies, we would like to press in that direction. It would be 
advantageous to spend entire nights observing lo continuously at 4.8 pm in 
order to obtain nonallased power spectra to 50 times lower frequencies. But 
such an undertaking %iould be at the expense of observing lo at other 
wavelengths and observing standard stars. With the data in hand, we are 
thwarted in doing a power spectrum by this problem of aliasing. One cannot 
numerically filter data that are incompletely sampled. In pressing toward the 
lowest frequencies that we shall discuss, we encounter numerous practical 
constraints that lu fact prevent complete sampling. Such constraints are (i) 
the necessity of looking at the same hemisphere of lo, (ii) lo can only be 
observed when it is above the horizon, (ill) Jupiter and clouds interfere some 
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times, and (iv) talcacopa achadullng coaaictaas do not hava unlinltad tlma to 
award. Thua av«n a program that la optimally daaignad to raduca allaaing will 
not do appreciably better than the data In hand tor the lowest frequencies at 
which we now know that lo fluctuatee. 

In the next section we will show that it 1s posulble to fit the available 
fluctuation data with a simple spectral dependence while using this spectral 
dependence to account for the aliasing Inherent In our data sampling. 


III. FLUCTUATIONS AT LOWER FREQUENCIES 

a. The Broadband Nature of the Fluctuation 

We are Impressed by two facts: %re are not able to find any repeatable 

maxima but rather find a trend toward lower frequencies In the power spectra, 
and the amplitude of flickering In the 2.2-m data (0.060 ± 0.005 fractionally) 
Is of the same order as the large night- to-night fluctuations seen In 4.8-pm 
photometry (Slnton^id. , 1983). Could the flickering and the nlght-to-nlght 
fluctuations be part of the same broadband phenomenon? Thus It Is Important 
to determine the fluctuations at an Intermediate time scale. To this t.?d. on 
9 and 10 June 1982 we acquired data In which we tried to make the 16 pairs of 
10-second Integrations (requiring 7.5 minutes %rlth the required 2-second 
settling time) every half-hour, along with a short 3.8- and 2.2-pm measurement 
and the same sequence on a standard star sandwiched Inbetween. We also 
searched our available data for sequences that could be used. Note that by 
forming the mean of each 7.5-minute segment, we average over rapid variations, 
and at the same time, the standard deviation determination for each night 
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•lialnatcs fluctuation* which ara vary alow coaparad to tha 4-5 hour* during 
which tha data la collactad* Tha praciaa diacuaaioo of tha affactlva 
frequancy bandwidth* la in a latar aactioo* 

All togathar 12 oighta of data ara availabla for thia analyaia. Those ara 
shown in Tabla II. As bafora, we hava to know tha aaacurad fluctuation of 
unvarying atars. This, detarmioad froa 34 sets of ■aasureaaots of stars on 
the saaa nights, is 0.035 fractionally. Tha oat ras fluctuation of lo is 
0.069 ± 0.013 of its total flux, which shows that lo varies on a tiae seal* 

that averages 40.0 ainutes between observations. The aaplitude is coaparable 
to that found \fith the fast flicker aaasureaaots, yat we have averaged each 
point over 7.5 ainutes of data. 

We can analyse our data for still lower frequencies. If, however, we 
coapare data froa adjacent nights, we are coaparing different sides of lo 
which have known differences. Hanes, we have coaparad nightly sveragas of 
data 48 hours apart, i.a., the saae side of lo. Thar* are six such pairs of 
nights listed lo Table III. Longer runs did not occur gcn*ra..iy because we 
request blocks of telescope tiae when lo is not close to Jupiter, and these 
blocks occur on a four-night cycle. We have assuaed an average tincertainty in 
the relative calibration of the paired nights of 5 percent. Thus the reduced 
res standard deviation of lo is 0.139 t 0.036, also expressed as a fraction to 
the total flux. As expected, lo shows a aarkad fluctuation on this tiae 
scale. We aay look at longer tiae scales and hence lower frequencies for 
f l^ictuations. We are, however ( constrained by the saapllng in the data in 
hand, and we cannot find longer repetitive Intervals in tha saapled data. The 
next longer asMnabl* tiae scale aayb* tha yearly averages. In Table IV of 
Sinton et (1983) there is a difference of 6 percent for the leading 
heaisphere between the two years. Considering that each yearly value is an 
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average of '>>80 neaaureiieata apread over ~20 alghta» arrora thla large are not 
expected. However, the aaae atar.dard atara have not bean ayatematlcally uaed 
and there la the poaalblllty of dlffarancaa avao anong groupa of the 
atandarda. We are preaently aeeklog to reaatahllah atandarda and reduce thla 
uncertainty. Nevertheleaa , It appaara that lo la variable from year to year. 
The program la continuing and by next year we will have doubled the time baae 
and will have tied together 45 atandard atara of aelected apectral typea to 
Improve the calibration. From 17 nlghta of obaervatlona obtained In 1982, the 
mean brlghtneea of lo for both the leading and trailing hemlapherea haa 
dropped by ~20 percent alnce 1980. Thus there appeara to be strong evidence 
of a year-to-year variation. 

b. Teat of Inverse Frequency Power Dependence 

Obviously, the fluctuation spectrum of lo Is very broad. Because of the 
aliasing, we also have very broad frequency response bands Inherent in the 
sampling techniques that tiere employed In Tables I-III. First, ve will make 
an estimate of the frequency dependence of the power spectrum. We will assume 
a spectral power dependence, and for each data set calculate how the aliasing 
will affect the response in a given frequency band. With the three different 
data seta we can then test If the assumed spectral dependence Is consistent 
with the response that was corrected for aliasing. 

The spectrum shown in Fig. 4 shows a definite tendency for the amplitude 
to increase with decreasing frequency. In reviewing the data acquired over 
various sampling Intervals In Tables I-III, we are Impressed by the fact that 
the standard deviation remains approximately 10 percent regardless of the 
sampling time. Despite our Inability to determine the power spectrum that Is 
free from the problems of aliasing for the longer period samplings, a little 
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reflectioo vlll show that the spectral aaplltude aust Increase for frequencies 
below those of the power spectral analysis already discussed. Let's suppose 
that lo had a white noise fluctuation, l.e., unlfora power at all frequencies. 
Then as we Increase the Integration tine, t, of a saaple we would expect the 
estimate of the standard deviation of the mean to decrease as The 

differences between samples aust also decrease as the precision of the 
estimates Improves. Yet It does not, even though the total Integration times 
Increases from ~20 seconds to thousands of seconds. In order for this to be 
true, the amount of the fluctuation of lo must Increase for the lower 
frequencies thac are being sampled as the Integration time and concomitantly 
the time between samples are Increased. 

From Just the observed constancy of the standard deviation. It Is 
possible to derive an estimate of the frequency spectrum of the fluctuation. 

As the total time for the array of measurements T Is Increased, the lowest 
frequency to .ruse a fluctuation between members of the array f decreases as 
1/T. If T always bears the same ratio to T, then the fluctuation for white 
noise «K>uld decrease as l/T^/^. Since In fact the observed fluctuation stays 
constant, we guess that the noise In a constant bandwidth must Increase as 
T^/2 or as f“l/2. 

It Is customary to express noise measurements lo terms of power. By 
power, the square of the amplitude of the fluctuation Is meant. Thus this 
simple estimate from the near constancy of the fluctuation gives a spectral 
dependence for the noise power within a unit bandwidth that varies Inversely 
as the frequency. It Is shown In Appendix A that the variance Is equal to the 
sum of the power over all of the frequencies that contribute to the 
fluctuation. As Just discussed, the lowest frequency Is 1/T. For data that 
are uniformly sampled In time one can obtain the power spectrum of the 
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fundanental and all harmonica • Paraeval'a relation (Appendix A) telle ua that 
the variance la equal to the total power at all frequenclea. 


n 



where la the power In the ^th frequency. If we make the aaaumptlon that 
the power spectrum « f~^ , then the summation can be written 
2 " -1 

- P, Z 1 , (9) 

* 1 - 1 

where we have made use of the equivalence of variance and total power and 
where Pj Is the power In the fundamental. It Is well-known that the series In 
(9) does not converge as n . However, the Integration of signal over 28^ 
In the case of Table I and 450^ In the case of Tables II-III effectively 
attenuates frequencies that lie above 1/56 Hz for Table I and 1/900 Hz for 
Tables II-III because for these frequencies the flux from lo was Integrated 
over more than 1/2 cycle. Thus ve can let n stop at an Integer Just greater 
than T/S6 for Table I and T/900 for Tables II-III. These Integers are 9, 16, 
and 192 respectively for Tables I-III, and the sums, S ■ r 1/1, In (9) are 
2.829, 3.381, and 5.837. The bandwidth, Af for each frequency Is T~1 . Thus 
the power per frequency bandwidth for each fundamental period Is given by 

Pi - o T/S . (10) 

These values are plotted as open circles in Fig. 4. A line with the slope of 
-1 Is shown In the logarithmic plot of this figure. 

What we have done, then. Is to assume for each set of data In Tables 
I-IIl that the noise power varies as f~l . Using that assumption and the 
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equlvaleoce of the variaoce to the total power, we have determined the 
amplitude of the lowest frequency for the data in each of the three tables. 
Comparison of these amplitudes in Fig. 4 to the line verifies the 
appropriateness of the f'^ assumption for the power spectrum. We have not 
sho%m that other spectral dependences cannot also work. A line with a slope 
as steep as -1.3 would also fit within the error bars of the extreme points. 
A slope of -1 is at the opposite extreme. 


IV. PHYSICAL MODEL FOR THE FREQUENCY DEPENDENCE 

The reciprocal dependence on frequency of the square of the fluctuations 
brings to mind several cases in which similar noise is developed in electronic 
apparatus. These Include many varieties of semiconductor devices and barium- 
oxide-coated cathodes of vacuum tubes. The theory of the latter has been 
elaborated (Schottky, 1926; Sproull, 1945; Macfarlane, 1947) and seems most 
analogous to the present case since electron micrographs of thermionic- 
emitting cathodes show that the emission occurs in evanascent patches that are 
called "hotspots." We can draw an analogy to the theory of the flicker effect 
developed by Macfarlane if we replace the thermionic current by the flux of 
thermal radiation. Note that in the electronics case, noise power is 
proportional to the square of the current fluctuations and that we find that 
the square of the thermal radiation fluctu.*itlons varies Inversely as the 
frequency. 

In the case of BaO cathodes, the theory supposes that clusters or patches 
of barium atoms diffuse to the surface, which Macfarlane termed "erupt," at a 
relatively small number of places. The presence of these atoms forms a 
polarized layer which locally reduces the work function for thermionic 
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emission. As electrons are emitted barium atoms are converted to Ions, so the 
polarization of the layer and the thermionic emission decrease ulth time. In 
drawing the analogy with thermal emission from lo's volcanoes » we consider 
that convection In cells In calderas will bring hot magma up from beneath. As 
the magma radiates Its heat away, the Initially high flux of thermal radiation 
will decay with time. Thus we replace diffusion of clusters of Ba atoms to 
the surface with diffusion of hot convection cells of a magma to the surface. 
We can write an equation for the rate of change of excess energy per cm^ , E, 
in the surface layer: 


dt 



( 11 ) 


where F is the flux of radiant energy from the surface, D Is the cellular 
diffusion coefficient, and h Is the distance over which the diffusion takes 
place. This equation In of the same form as Eq. (1) of Macfarlane (1947). 

The other basic equation used by Macfarlane Is Langmuir's equation which 
relates the current to the excess of atoms on the surface. The equivalent 
equation relating F to E can be found from the relation of the surface energy 
content to the specific heat, C, and temperature, T, E > hCT; and from the 
Stefan-Boltzmann law, F > a T*' . Combining these gives 


F 


^ 4 4 • 

h C 


( 12 ) 


This does not have the same functional form as Langmuir's equation, tfhich Is 
exponential In the number of atoms. The author has solved the differential 
Eq. (11) when combined with (12) exactly, while Macfarlane found only an 
approximate solution to his equation. The two solutions, however, have 
approximately the same time dependence. 
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One additional assumption made by Macfarlane was that the number of 
excess atoms had a Gaussian distribution with an e-£oldlng width of N, which 
would be analogous to a similar distribution In the Initial temperature of the 
convection cells on lo. We will not repeat the lengthy mathematical 
devel'^pment here* However* he shows that the dependence of noise power on 


current* J* and frequency follows the relationship 


< A J > 


1 


1 + X 


(13) 


where 1 < x < 2 for « > N > 0. The empirical value of x for cathodes Is 
1.25* In the previous section we found an empirical value of x for lo which 
lies between 1.0 and 1*3. 

The similarity of the formulation of cathode flickering to what Is known 
about the volcanlsm of lo and the spectrum of the fluctuation that Is 
reported here make the development of the flicker theory for lo look very 
attractive. The author Intends to make this application and Investigate the 
implications of this theory for lo volcanlsm In a future paper* 


V. SPECULATION ON LONG-TERM FLUCTUATIONS 


Evidence was previously given that lo fluctuates from year to year on a 
global scale at 4.8 ^m* In this section we will speculate on the consequences 
of variability at even longer time scales* say 100 years or more. Eclipse 
measurements have shown that •^40 percent of Io*s flux at 4*8 pm Is thermal 
emission from volcanoes (Slnton et al. 1980). There Is already evidence that 
the flux varies by ~20 percent from year to year. Since these are similar* It 
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teems reatoaable to suppose that excursions may reach 40 percent from the 
current value and that not only can the activity double, but It can also cease 
completely and In fact may do so for very long times* The times might be long 
enough that the present vents will freeze and thin places In the crust 
thicken. Thus lo may go from Its present volcanically active state to an 
inactive one. Slnton (1982) discusses this possibility as an on-off 
thermostatic mechanism. The fluctuations discussed here and the apparently 
open-ended spectrum derived thus far constitute the first Indication of how 
this might come about* Furthermore, the measurements of total volcanic heat 
flow (Matson et al. , 1981; Slnton, 1981; Morrison and Telesco, 1980; Pearl and 
Slnton, 1982) when compared to estimates of tidal dissipation (Yoder, 1979; 
Cassen et al. , 1982; Greenberg, 1982; Yoder and Peale, 1981) suggest that lo 
is currently losing more heat than It la gaining. Hence In the long term, 
activity may cease temporarily because lo may not be able to maintain the 
current level. 


VI. FLICKERING TESTS AT 10 ;im 

It seems obvious that the flickering at 4.8 ^m results from ongoing 
volcanic activity on lo. In an earlier paper (Slnton, 1982) we presented a 
steady-state model to explain the spectrum of lo observed from 2*2 to 30 ^m 
during an eclipse. In the model we ascribed the thermal emission at the short 
wavelengths (2. 2-4. 8 ^m) to current eruptions, while the radiation at longer 
wavelengths was ascribed to the cooling of "old" flows. To test this Idea, we 
reasoned that old, solidified flows would not flicker. Fig. 4 shows rapid 
time series measurements of lo at 10 outside of eclipse. Flickering with 
an amplitude of 1 percent Is not detectable at 4o in these data. 
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However, at this point the negative result for flickering should not be 
cciisldered as definitive because lo was not measurably flickering at 4.8 
either. More observatlona are deemed necessary for nonfllckerlng to be 
considered as evidence for the flow model. 

Eclipse observations at 10 pm have thus far given quite consistent 
results. However, If the low flux at 4.8 pm observed In 1982 Is maintained, 
we should expect a decrease In the 10 pm flux since, according to the flow 
model, the current activity as viewed at 4.8 pm produces the flows that are 
subsequently seen at 10 pm. Thus In future work we should cross-correlate the 
fluctuations observed at 4.8 with any variability measured at 10 pm as a 
further test of the flow model. 


VII. CONCLUSIONS 


It has been shown that the 4.8-pm flux from lo Is variable at any time 
scale at which lo has been observed. The spectrum of this fluctuation Is such 
that the mean square fluctuation In a unit frequency Interval varies Inversely 
as the frequency within the errors of the determination over a frequency range 
of 6 pHz to 17 mHz. This type of fluctuation or "noise" Is similar to many 
sources of noise found In electronic equipment. A theory that has been 
developed for BaO cathodes appears to have many similarities to the volcanlsm 
on lo. Including the appearance of evanescent "hot spots" on the surface and 
the diffusion of clusters or cells which transport that which Is to be emitted 
to the surface. In the case of cathodes, electrons stored In Ba atoms are 
transported; In the case of lo, the heat stored In magma Is transported. The 
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nuithematical traataant la IdanClcal la tha dlffualoo-coaducclon aqtatloo, but 
diffara In tha relation batwean atoraga and emlaalon. 

We speculate that the fluctuation of tha volcanic activity on lo may 
extend to very long time parloda» and there la evidence that It extenda at 
least to Intervals of several years. Such very low frequency variation may 
offer A mechanism whereby the volcanlsm can be Intermittent as the author has 
suggested in a previous paper. 

Lastly, It was found that on one night In which tests were made, lo 
did not flicker at 10 pm. This tentative result la In agreement with 
predictions of the author's "flow model," but further tests are required. 
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APPENDIX A 

RELATIONSHIP BETWEEN THE TOTAL SPECTRAL POWER AND THE VARIANCE 

We have not been able to find thla relatloaahlp lo any popular textbook 
probably becauae moat textbooks treat Fourier series separately from 
statistics. Parseval's relation for Fourier series gives 

ci2 - |/T [ f(t)]2dt 
^ 1-1 ° 

where ^ f(t) dt 

and Cl} - ai^ ^ b^2 • ^ 2 f(t) cos t ) dt 

+ [ I /j *(t) ‘lo (-^ t ) -It • 

since cj^ is the amplitude of the sinusoidal wave of harmonic 1* the power 
in the harmonic Is 


2 



where the factor of 1/2 arises from taking the average power over a sinusoidal 
wave. The constant term, Sq/2, Is the average value of f(t) over the Interval 
0 < t < T. If the function is sampled at n equally spaced points over the 


interval we may write 
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^ '(f)- 1“ i ? ,'3 

n ♦ • J«1 n ♦ • J«l 

Approximating the Integral by a aua using the trapezoidal rule, 

i/J[f(t)]*dt- lU I ? 11. I I 

0 n • J"1 n ♦ • J"1 


me have 


2 


Thus we have: 


n 

Z P, - 11m 
1-1 n • 






2 


In the limit as n • the left side la the total power In all the harmonics 
while the right side Is the variance. It may be seen that In an electrical 
circuit the left side expresses the total power In all the ac harmonics while 
the right side expresses the total power as proportional to the mean square 
voltage or current variations. Conservation of energy requires that the 
equality hold. 
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TABLE I 

Fluctuation of Signal Ralativa to Total Signal at 28-Sacond Intarvala 


Data 

Hi 

Star 

»b 

lo 

oa 

Data 

lo 

a 

Star 

Ob 

o» 

is fak sr~ 




25 Fab 81 





0.026 

0.065 

0.060 ^ 


0.062 

0.050 

0.037 


0.070 

0.065 

0.026 


0.057 

0.050 

0.027 


0.085« 

0.065 

0.055 


0.067* 

0.050 

0.045 


0.121* 

0.065 

0.102 


0.055 

0.053 

0.015 


0.033 

0.065 

0.038 b 


0.087* 

0.053 

0.069 


0.058 

0.065 

0.029 


0.182* 

0.123 

0.134 


0.131* 

0.065 

0.114 


0.168* 

0.069 

0.153 


0.064 

0.065 

0.011 ^ 


0.193* 

0.120 

0.151 


0.067 

0.065 

0.016 


0.091 

0.069 

0.059 


0.059 

0.065 

0.027 *> 


0.156 

0.123 

0.096 


0.062 

0.065 

0.020 


0.186* 

0.123 

0.140 

24 Fab 81 

0.114* 

0.055 

0.100 


0.163 

0.123 

0.107 


0.134* 

0.048 

0.125 

19 Mar 81 

0.107 

0.064 

0.082 


0.035 

0.046 

0.030 b 


0.067 

0.069 

0.016 *> 


0.095 

0.050 

0.083 


0.089 

0.069 

0.056 


0.037 

0.046 

0.027 b 


0.069 

0.069 

0.000 


0.046 

0.046 

0.000 


0.086 

0.069 

0.051 


0.047 

0.046 

0.017 


0.115* 

0.069 

0.092 


0.070 

0.046 

0.053 


0.128 

0.080 

0.100 

23 Feb 81 

0.044 

0.048 

0.019 *> 


0.164* 

0.069 

0.149 


0.053 

0.053 

0.000 

2 Jan 82 

0.138* 

0.055 

0.127 


0.061 

0.055 

0.040 


0.148* 

0.058 

0.136 

23 Feb 81 

0.048* 

0.053 

0.022 »> 

4 Jan 82 

0.080* 

0.051 

0.062 
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Table I (continued) 


Date 

lo 

a 

Star 

lo 

oe 


Date 

lo 

a 

Star 

Ob 

lo 

Os 

4 Jan 82 




4 Feb 82 

0*050 

0*034 

0*037 


0*087« 

0*051 

0.077 



0.035 

0*034 

0*008 


0*092« 

0*051 

0*077 



0*030 

0*034 

0*016 ^ 


\ 

0*097* 

0*051 

0*083 



0*051 

0*034 

C.038 


0*097* 

0*051 

0*083 



0*048 

0*034 

0*034 


0*091* 

0*051 

0*075 



0*028 

0*034 

0*019 b 


0*090 

0*051 

0*074 



0*030 

0.034 

0*016 b 

5 Jan 82 

0.063 

0*051 

0*037 



0*042 

0*034 

0*025 


0*065 

0*051 

0*040 



0*028 

0*034 

0*019 b 


0*055 

0*051 

0*021 



0*038 

0*042 

0*018 b 


0*059 

0*051 

0*030 



0*041 

0*042 

0*009 b 


0*064 

0*051 

0*039 


5 Feb 82 

0*048 

0*042 

0*023 


0*049 

0*051 

0*014 b 



0*033 

0*042 

0.026 b 


0*093* 

0*051 

0*078 



0*028 

0*042 

0*031 b 


0*055 

0*051 

0*021 



0*042 

0*042 

0*000 


0*065 

0*051 

0*040 



0*050 

0*042 

0.027 

11 Jan 82 

0*069 

0*049 

0*049 



0*040 

0*042 

0.013 b 


0*056 

0*049 

0*027 



0*048 

0*042 

0*023 


0*051 

0*049 

0*014 



0*039 

0*042 

0*016 b 


0*093* 

0*049 

0*079 



0*042 

0*042 

0*000 


0*055 

0*049 

0*025 



0*033 

0*042 

0*026 b 




/ 



0*049 

0*042 

0.025 


Root-mean-square lo Cast flicker, N/S ■ 0«060 ± 0*005 • Flicker in 
magnitudes * 0*062 ± 0*005* 

^ These data sets were used for the computation of the power spectrum* 

^ These values are Imaginary and are Included lo the cniculatloo of the 


root-meea-square flicker 
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TABLE II 

Fluctuation of Signal Belatlve to Total Signal 
at Noalnal aO-Mlnute Intarvala 


Date 

Number 

Average 

Spacing 

(min) 

a 

Os 

23 Feb 81 

9 

37.5 

0.067 

0.056 

24 

6 

36.0 

0.151 

0.146 

25 

9 

39.5 

0.049 

0.035 

18 Mar 81 

8 

58.6 

0.137 

0.131 

19 81 

7 

62.0 

0.089 

0.081 

4 Jan 82 

5 

33.5 

0.024 

0.025 « 

5 82 

5 

46.2 

0.039 

0.019 

11 82 

3 

31.5 

0.058 

0.045 

4 Feb 82 

7 

33.2 

0.014 

0.032 a 

5 82 

6 

32.2 

0.062 

0.056 

9 Jun 82 

9 

32.0 

0.070 

0.060 

10 82 

5 

38.5 

0.039 

0.019 


Measured fluctuations of stars on the same nights from 34 sets of 


2 to 6 sKasurenents 0*035 mag. 

‘These numbers are Imaginary and are Included In the calculation of the 
root-nean-aquare fluctuation. 

Root-mean-square lo fluctuation 0.069 i 0.013. 
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TABLE III 

Observations of the Same Hemisphere Separated by 48 Hours 


Hemisphere 

Date 

Duration Number 

Mean 

8 • d • 

Corrected 

for 


UT 

min. 

2.5 log p® 
mag. 

mag. 

Measure- 

ment 

error 

mag. 


Leading 

23 

Feb 

81 

308 

13 

-0.096] 









J 

0.172 

0.165 


25 



329 

12 

0.147] 




2 

Jan 

82 

19 

2 

-0.250] 









1 

0.089 

0.074 


4 


82 

134 

5 

-0.124] 



Trailing 

13 

Dec 

79 

301 

7 

0.250] 

0.271 

0.266 


15 



178 

6 

-0.133] 




20 



249 

5 

0.110] 




22 



346 

13 

-0.120' 

< 

0.125 

0.115 


24 



— 

1 

0.08 ; 



Perl jovlan 

22 

Apr 

81 

44 

3 

-0.020] 









j 

0.071 

0.050 


24 



63 

6 

-0.120] 



Apo Jovian 

21 

Apr 

81 

98 

3 

-0.203] 










0.057 

0.027 


23 


81 

286 

12 

-0.123] 




^The quantity p Is the quasl-geometrlc albedo as defined In Slnton et . al 
(1983). 

Net rms variation of lo ■ 0.141 ± 0.038 mag. 

Variation reduced to fraction of signal ■ 0.139 ± 0.036. 
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FIGURE CAPTIONS 

Fig. 1. The aolse-to-slgnal ratloe, N/S, for stars (points) and lo (open 
circles) determined during various observing runs at 4.8 pm. The curves are 
weighted linear regression fits to the stellar data In the form of Eq. (8). 

The observing runs for several panels are (a) 22-25 Feb 81 » (b) 18-19 Mar 
81 » (c) 2-5 Jan 82, (d) 10-12 Jan 82. The error bars for all points except 
some of those of panel (a) are (1/30)^^^ of the value. For lo measurements 
that lie above the dotted curve, there Is 95 percent confidence that lo was 
flicker* 

Fig. 2. The nolse-to-slgnal ratio for the observing run 3-5 Feb 82 for 
stars (open circles) and lo (11 points In the speckled box). During this 
run lo flickered very little. 

Fig. 3. Rapid changes In lo's flux at 4.8 pm, 21 Apr 81 (upper) two 
traces. The lower traces of stars show the IRTF bolometer system noise. All 
measurements were made %d.th the same gain and the stars have nearly the same 
brightness as lo. The breaks In the traces were to check centering. 

Fig. 4. The power spectrum of Io*s fluctuations. The solid points are 
from the mean power spectrum derived from the selected 21 time series. The 
plotted points have been reduced so that selection effect Is eliminated, and 
their average power Is the same as for the entire 91 rapid time series. The 
open circles are spectral points derived by the method described In the text 
from sequences every <^40 minutes (middle point), every 48 hours (left point), 
and the 91 rapid time series (right point). The line has a slope of -1. 

Fig. 5. Time series of measurements of lo's flux at 10 pm, 24 Apr 81 
(upper three traces). The star observed In the lower trace was about three 
times brighter than lo, but the data are shown at the same gain. The drift 
evident In the star's trace Is caused by a slight tracking error. 
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